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ABSTRACT 

Energetic electrons of up to tens of MeV are created during explosive phenomena in the solar 
corona. While many theoretical models consider magnetic reconnection as a possible way of generating 
energetic electrons, the precise roles of magnetic reconnection during acceleration and heating of 
electrons still remain unclear. Here we show from 2D particle-in-cell simulations that coalescence 
of magnetic islands that naturally form as a consequence of tearing mode instability and associated 
magnetic reconnection leads to efficient energization of electrons. The key process is the secondary 
magnetic reconnection at the merging points, or the 'anti- reconnection', which is, in a sense, driven by 
the converging outflows from the initial magnetic reconnection regions. By following the trajectories 
of the most energetic electrons, we found a variety of different acceleration mechanisms but the 
energization at the anti-reconnection is found to be the most important process. We discuss possible 
applications to the energetic electrons observed in the solar flares. We anticipate our results to be 
a starting point for more sophisticated models of particle acceleration during the explosive energy 
release phenomena. 

Subject headings: acceleration of particles — Sun: flares — Sun: X-rays, gamma rays 



1. INTRODUCTION 

A solar flare is an explosive energy release phenomena 
on the sun and a large fraction of the released energy ap- 
pears to go to high ener gy, often non-th ermal, particles 
both ions and electrons (|Lin et al1l2003L and references 
therein) . The particle energy reaches tens of GeV for ions 
and tens of MeV for electrons. The mechanism of pro- 
ducing such energetic particles is much less understood 
compared to the energy release mechanism. Because 
magnetic reconnection is a possible mechanism of the en- 
ergy release process, particle acceleration may also occur 
through magnetic reconnection, although difficul ties re- 
maiii when trying to interpret observatio ns (Miller et al.l 
Il997t iKrucker et all 120081: iBenz 1 120081 and references 
therein) . In order to explain obs ervations, many the oret- 
ical ideas have been proposed ()Aschwanden I l2002l and 
references therein) . While some theories utilize fast / slow 
mode shocks as well as electromagnetic waves of various 
scales - sometimes in a stochastic manner - many theo- 
ries consider magnetic reconnection. In this paper, we 
also assume a priori that magnetic reconnection plays a 
role for particle acceleration and explore possible mecha- 
nism(s) of particle acceleration in association with mag- 
netic reconnection. 

In general, a test particle approach had been used 
to explore particle acceleration by magnetic reconnec- 
tion. While some studies solved particle motion ana- 
lytically (e.g. iLitvinenkol 11996') . a majority of studies 
performed test-particle simulations under model electro- 
magnetic fields (e.g. [Kli cni 1994; Hannah & Fleichex 
[200ft lOnofri. IslikeTfcVlahos I |2006|) or time varying 
fields generated by MHD simulations (iSato et al. II1982I : 
iScholer and Jamitzk"v1ll987l:lAmbrosiano et al.lll988[ ). A 
self-consistent, particle-in-cell (PIC) simulations are now 
widely used to study the detailed process of electron 
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energization , (e.g. iHoshino et al~l 120011: iHoshino I 200 



Drake et al.' 2005; Pritchett ' 2006; 'Drake et al. 200€ 
Karlic kv & Barta 2007; Pritchett 2008; Wan et alL 
2008t rShinohara et al.l 1200911 . PIC simulations have an 



advantage of being able to resolve the inner structure 
of the X-line, the so-called diffusion region. It is a sci- 
entific challenge to understand particle acceleration by 
magnetic reconnection that involves multi-scale. 

For convenience, we categorize various theories of par- 
ticle acceleration associated with magnetic reconnection 
into two different groups: X-type acceleration and 0-type 
acceleration. The X-type acceleration takes place at and 
around the X-line of magnetic reconnection or the dif- 
fusion region. Particles are unmagnetized at the X-lines 
and can be direct ly accelerated by the electric field (e.g. 
iSato et al. II1982I ). In the immediate downstream of the 
X-line are the regions with magnetic field gradient where 



rent sheet (e.g. 


Scholer and Jamitzkv 1119871 iKliem II1994I: 


IHoshino et al. 


I2001D. More recentlv. it was found that 



the in-plane, polarization electric field in the diffusion 
region generated by the charge separat ion between ion s 
and electrons plays an important role ([Hoshino 1 120051) . 
The force by the polarization electric field can be bal- 
anced by the Lorentz force so that electrons are acceler- 
ated efficiently by the reconnection electric field while be- 
ing trapped within the current sheet boundary. Because 
of the trapping effect, the process was named 'surfing' 
mechanism. 

The 0-type acceleration takes advantage of the closed 
geometry of field lines in a magnetic island. In many 
cases, a magnetic island is bounded by two X-lines at 
each end. Therefore, if particles are trapped in a mag- 
netic island, they can continue gaining energy by repet- 
itive crossings of the gradient region, although the re- 
connection el ectric field is relatively weak inside mag- 
netic islands (jStern 1 119791 IScholer and Jamitzkvl 119871 : 
IKliem |[l994[ ). In order to compensate for this weak elec- 
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Reoonnectkjn 
by the tearing instability 





Anti-Reconnedion (reversed E) 
by the island coalescence 

Fig. 1. — Schematic illustration of multi-island coalescence. The 
thin cross marks indicate the X-lines of normal magnetic reconnec- 
tion. The thick cross mark indicates the X-line of the anti-magnet- 
reconnection generated by the coalescence. The arrows indicate 
the flow directions. 



trie field, a recent model has been developed that takes 
into account the dy namical, contracting motion of islands 
(jPrake et al.l 120061 ). The time-dependent model is anal- 
ogous to the energy increase of a ball reflecting between 
two converging walls, namely the first order 'Fermi' pro- 
cess. 

Despite the intensive research of magnetic reconnec- 
tion, it still remains unknown which of the two differ- 
ent type of acceleration is important for producing en- 
ergetic electrons. In this respect, multi-island coales- 
cence has drawn considerable attentions because it po- 
tentially contains both X-Type and 0-Type mechanisms. 
A coalescence is the process of merging of two mag- 
netic islands and has been studied by both theories and 
simulations (Finn & Kaw' '1977^ 'Pritchett & Wu ' '197St 
Biskam p fc Welte r 1980; Tajima ct al. 1987; Pritchett! 
2007t (Wan et al. I [20081. Figure 1 shows a schematic il- 
lustration. In general, multiple number of X-lines and lo- 
calized currents are generated from the tearing-mode in- 
stability (Figure[T^). If we assume an X-line at the center 
was relatively weak compared to the other two X-lines, 
the two localized currents bounding the central X-line 
will eventually be attracted to each other by the Lorentz 
force (Figure [T})) . Each current is represented by mag- 
netic islands and at the merging point of the two islands, 
a secondary mag netic reconnection or 'anti-reconnection' 
(jPritchett Ii2"00l occurs (Figure[TJ;). The direction of the 
electric field of the anti-reconnection is reversed from the 
direction of the electric field of the primary reconnection. 
Finally, the two magnetic islands become one large mag- 
netic island. 

While the first PIC simulation of particle acceleration 
during multi-island coalescence was performed more than 



decades ago (iTaiima et al. II1987D. a d etailed study came 
out quite recently. iPritchett I ( 20081) showed that elec- 
trons are energized as the number of magnetic islands is 
reduced by coalescence. An important conclusion of the 
study is that the reversed electric field decelerates elec- 
trons near the anti-reconnection site. It was suggested 
that the main energization occurs through the 0-Type 
mechanism. 

In this Paper, we extend the work by IPritchett! ()2008f ) 
by performing 2D PIC simulations of multi-island coales- 
cence with no guide field. The key and rather surprising 
finding of our simulations is that the anti-reconnection 
plays an important ro l e in a ccelerating electrons. Very 
recently, iTanaka et al.l (|20100 also reported intense elec- 
tron energization by the anti-reconnection, but our work 
provides new insights from different perspectives on this 
issue because we fully analyzed the trajectories of ac- 
celerating electrons and clarified the importance of the 
anti-reconnection with respect to the other mechanisms. 

The outline of the paper is as follows. In section 2, 
we describe the setup of the simulation runs and the 
overview of results. Section 3 presents analysis of elec- 
tron energy spectra. Section 4 presents the trajectory 
analysis of energetic electrons and describe the details of 
their energization processes. In section 5, we summarize 
the results and discuss applications to the observations 
of the solar flares. Finally, section 6 contains the conclu- 
sion. 

2. SIMULATION SETUP AND OVERVIEW 

We utilize a two and half dimensional, fully rel ativis- 
tic PIC code (lHoshinolll987l:IShinohara et al.l[200l . The 
initial condition consists of two Harris current sheets. 
The anti-parallel magnetic fleld is given by By/Bo = 
tanh((a; — xii)/D) — tanh((a; — xl)/D) — 1 where Bq is 
the magnetic fleld at the inflow boundary, D is the half- 
thickness of the current sheet and and Ly are the 
domain size in x and y direction, respectively. xl—Lx/4: 
and xii=3Lx/4: are the x-positions of the left and right 
current sheet, respectively. Periodic boundaries are used 
in both directions. The ion inertial length di is re- 
solved by 25 cells. D=0.5d, and = Ly=l02Ad^. The 
inflow, background plasma has the uniform density of 
A^BO=0.2A^o where A^o is the density at the current sheet 
center. The ion to electron temperature ratio is set to 
be Ti/Te=5 for the current sheet and Ti/Te=l for the 
background. The background to current sheet tempera- 
ture ratio Tiik/Tcs—O-l- The frequency ratio u)po/^cc=S 
where ujpc and Qcq are the electron plasma frequency and 
the electron cyclotron frequency, respectively. The ion 
to electron mass ratio TOi/me=25 and the light speed 
c is 15V A wh ere Va is the Alfven speed defined as 
Bq / V47rAomi. We used average of 64 particles in each 
cell. 297 particles per cell represents the unit density. 
No magnetic field perturbation is imposed at the begin- 
ning so that the system evolves from the tearing mode 
instability due to particle noise. 

The overall evolution of the sys tem is v ery much the 
same as those of previous reports (jPritchet t 2008). The 
evolution of our simulation run is summarized in Fig- 
ure [5] and [3] The linear growth of the tearing instability 
leads to the generation of at least seven X-lines by the 
time i7cit=43 (Figure [2^). In the non- linear stage, the 
number of X-line is reduced due to coalescence (Figure 
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Fig. 2. — Out-of-planc component of the electron current density obtained at (a) Slcit = 42.7, (b) licit = 70.0, (c) f^cit = 102.0, (d) 
Qcit = 114.7. The fihed circles are the positions of the energetic electrons with energies (c) e >1.2moc2 and (d) e >1.4mcc2. The arrows 
indicate island merging points. 



[2)D-d). At each merging point, anti-reconnection occurs 
and the out-of-plane component of the electron current 
density 3e,z enhances at, for example, (x, y)=(25,84) and 
(77,40) in Figure [2j;. Small magnetic islands by the sec- 
ondary tearing instability are also observed at, for ex- 
ample, (x, i/)=(26,31) in Figure but they will merge 
with a larger island. By the end of the simulation run, 
one large island remains in each current sheet. 

Note that if two current sheets are sufficiently close, the 
tearing mode and magnetic reconnection drive each other 
(|Pritchettl[l980l ). In the present case, the two current 
sheets evolve independently until ricit'~-"75. In the later 
stages of the simulation (r2cit>105), the configuration of 
one layer becomes almost anti-symmetric to the other 
layer. However, by varying the separation between the 
two layers (or L^), we verified that results presented in 
this paper are insensitive to the domain size. 

Figure[5^ shows the time profile of particle and field en- 
ergies integrated over the entire simulation domain. Co- 
alescence proceeds in two distinct stages: 40< fldt <65 
and 90< ildt J;115. During these intervals, the mag- 



netic field energy is more rapidly converted to particle 
energy. At the end of the run, ions carry 67% of the 
released magnetic field energy whereas electrons carry 
33%. About 34% of the initial magnetic field energy has 
been released. Figure ^jp shows the time profiles of the 
number of electrons of certain energies. It is apparent 
that higher energy flux increases rapidly in association 
with the two growth stages mentioned above. 

In order to explore origins of the energetic electrons, we 
plotted locations of the most energetic electrons in Figure 
[2J:,d. Figure corresponds to the time of rapid increase 
of the energetic electron flux. At this time, electrons are 
distributed near the merging point suggesting that these 
electrons are created by the anti-reconnection. Figure 
corresponds to the time of the energetic electron peak 
flux. By this time, coalescence is almost complete and 
energetic electrons are uniformly distributed inside the 
merged island but surrounding the core region. It is also 
important to note that energetic electrons also exist in 
the secondary magnetic island. 

3. ENERGY SPECTRUM 
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Fig. 3. — (a) Time histories of particle and field energies inte- 
grated over the entire simulation domain. The electric field energy 
profile has been multiplied by 1000. (b) Stack plots of the time 
histories of the number of electrons counted in each energy bin Si 
where Si = (0. 0125+0. 025i)mec2, i=0,l,2..., from top to bottom. 
The bin size A£i=±0.0125mcC^. Some of the bin energies are an- 
notated in the panel. The gray arrows in between the two panels 
indicate the time of the snapshots of Figure [2] 

Figure |4^ shows the time variation of the electron en- 
ergy spectrum integrated over the entire simulation do- 
main. Electron heating starts during the earlier phase 
of the initial growth phase (r2cit=43). By the end of 
the initial growth phase (ricit=72), the spectrum con- 
sists of at least two components: a cold thermal compo- 
nent representing the initial plasma condition and a hot 
thermal component representing the heated/accelerated 
plasma. This spectral form does not change very much 
there after but the heating/acceleration continues con- 
stantly (ficit=102, 115). Figure \^ shows the energy 
spectra integrated over the rectangular boxes in Figure 
[2j;. These are obtained when the highest energy electrons 
appeared around the anti-reconnection site (r2cit=102). 
It is shown that the hot, energized electrons already exist 
in the vicinity of the X-line (A) but the magnetic field 
pile- up region produces more energetic electrons (B). The 
most energetic electrons, however, are created in the anti- 
reconnection region (C). FigurelH; shows the energy spec- 
tra integrated over the rectangular boxes in Figure 
These are obtained when the flux of energetic electrons 
was at its peak (rJcit=115). It is evident that the most 
energetic electrons exist in the merged island (D), al- 
though the secondary island alone can also produce a 
hot thermal component (E). The above results suggest 
that the most energetic electrons are first accelerated by 
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Fig. 4. — (a) Time variation of the energy spectrum integrated 
over the entire simulation domain, (b) Energy spectra obtained at 
f2t.it=102 in the area A, B and C of Figure[2};. (c) Energy spectra 
obtained at f2cit=115 in the area D and E of Figuro[2jl. 

the anti-reconnection and then further energized within 
the merged islands. 

Figure [5] shows exactly the same spectrum shown in 
Figure |4^ obtained at ficit^HS, but it is now fitted by 
a best fit model: 

F{e) = Fo + A, exp(^-^^ 

+ ^2exp(-^)+A3exp(^-^) (1) 
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Fig. 5. — Electron energy spectrum obtained at f2cit=115. The 
gray curve shows the best fit model shown by the equation JTJ. 
The dashed, dotted, and dash-dotted curves show the first, sec- 
ond, and third thermal component of the best fit model. For the 
fitting procedure, we used the standard deviation of particle num- 
ber count in each energy channel to evaluate the weight. The 
resultant numbers are Fq = 1.12±0.15, Ai = 1.52x 10**±1.84x lO"*, 
Ti = 7.3xl0-^±1.3xl0-^moc2, A2 = 2.11xl0''±6.96xl03, 
T2 = 0.04±1.3xl0-5mcc2, A3 = 9.90xl05±2.75xl0=', T3 = 
0.11±8.4xl0-'^mcc2, and x^/d-o.f.=2.7x lO'^/SSO. For the hori- 
zontal axis, the energy e is normalized by the rest mass energy 
meC^ (bottom axis) and the initial background temperature £th 
(top axis). 

The first exponential component (Ti ^7.3x 10~^mcC^, 
dashed curve) corresponds to the initial electron 
distributions. The second exponential component 
(T2 ~4.2x 10"'^mcC^, dotted curve) represents the heated 
electrons by the multi-island coalescence and associated 
magnetic reconnection. The third exponential compo- 
nent (T3 ^O.llmcC^, dash-dotted curve) may be viewed 
as either additional thermal component or a non-thermal 
component that was not able to extend to higher energy 
possibly because of the limited size of the simulation. For 
a reference, we indicated a slope of 3.5 by the solid line. 

4. PARTICLE TRAJECTORY 

Knowing that multi-island coalescence produces a very 
hot, possibly non-thermal component, we focus on the 
highest energy electrons. In order to clarify the ex- 
act process of electron energization, we traced backward 
in time the trajectories of the electrons with energies 
e >1.4mcC^ plotted in Figure [2jl, that is, in total, 198 
electrons. Among these particles, we have found a variety 
of different energization processes and organized them by 
the type of acceleration mechanism. 

4.1. X-type Acceleration 

Figure m shows the first example of the electron trajec- 
tories (hereafter Particle #1). We divided the trajectory 
into four segments and plotted over the out-of-plane com- 
ponent of the electron current density. The whole profiles 
of energy and the displacement Az are shown in the right 
two panels. The simulation itself is two-dimensional in x 



and y but the displacement Az can be calculated by in- 
tegrating Vz over time. The background images are just 
snapshots taken during each segment so care should be 
taken when comparing the trajectories with the images. 
The time of each image is chosen so that it best reflects 
the characteristics of each trajectory. 

Particle #1 is first energized at one of the few X-lines 
formed within the first coalescence stage (Figure |6^). 
Scatterings take place in the immediate downstream of 
the X-line, but the particle soon starts to travel along 
the rim of the island associated with the X-line and is 
further scattered when it reaches near the merging point 
of two islands approaching to each other (Figure ^jp). 
While wandering around the merging point, the coales- 
cence progresses to form a region of localized current 
in which the particle is then confined to experience the 
second rapid energization (Figure |6l;) . Note that this lo- 
calized current corresponds to the anti-reconnection in 
which electric field direction is reversed. The energiza- 
tion persists even after the particle is ejected out of the 
central anti-reconnection region and the particle is pitch 
angle scattered afterwards (Figure |6ji). The first and 
second rapid energization occur in a very localized re- 
gion (Ay ~ 0.1 di) associated with the reconnection and 
anti-reconnection, respectively (Figure IBJ3) . What is im- 
portant here is that the particle moves toward the pos- 
itive z-direction up to Az ~50di, but by the reversed 
electric field, it turns its direction toward the negative 
z-direction. As a result, it passes through the original 
position and reaches Az ~ -lOdi by r2cit=104.0. 

Figure [7] shows the details of the behavior of Particle 
#1 during the first rapid energization. The trajectory is 
further divided into three time periods according to its 
characteristics. The time profiles of the particle energy 
as well as the electric fields felt by the particle are also 
displayed. 

Figure [7^ illustrates the period during which the par- 
ticle approaches toward the X-line and starts to be ener- 
gized. During this phase, we can identify the polarized 
electric field mainly parallel to the reconnection plane 
(i.e. X — Y plane) and the particle is energized as it 
passes through this region. The polarized electric field 
felt by the particle reaches ~ VaBq (Figure [7^). How- 
ever, the corresponding energy increase is very small, 
Ae ~0.1meC^ (Figure [fjl) . In order to better visualize 
the energy increase, we multiplied the energy profile by 
10, as shown by the dashed curve. 

Figure [Tb illustrates the 'X-line Phase' during which 
the particle stays at the very center of the diffusion region 
or the X-line and is rapidly energized. Within a relatively 
short period of time r2ciAt'-^2, the energy increases nearly 
an order of magnitude. Note also that this energization 
takes place in the so-called inner diffusion region. Recent 
simulations revealed that the out-of-plane current is lo- 
calized only at the center of the diffusion region ('inner 
diffusion region') so that the reconnection rate remains 
fast while electro ns form a high- velocity jet in the 'outer 
diffusion r egion ' (Daughton et al."2006; Fuiimoto 20061 
Shav et al.„ 2007; Karimabadi_et al. 2007; Phan ct al| 
2007D . Electrons are unmagnetized in these regions so 
that a snapshot of the z-component of E-|-Ve x B makes 
clear in which region the electron can be energized. The 
inner diffusion region (colored blue) is at the X-line and 
the outer diffusion region (colored red) extends toward 
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Fig. 6.— The trajectory of particle #1 during (a) 62. 0< Q^it <75.3, (b) 75. 3< fl^it <87.4, (c) 87.4< fl^it <93.8, (d) 93.8< Q^it <104.0. 
Also displayed are (e) the energy and (f) the displacement AZ during 0< fi^it <114.7 as functions of Y. The background images are the 
electron current density at (a) r2t.it=70.0, (b) ncit=84.0, (a) ncit=92.7, (a) ncit=104.0. 



the downstream. It is evident that the electron is rapidly 
energized within the 'inner diffusion region'. A close look 
at the particle trajectory reveals a meandering motion in 
this region (Figure [7] inset). 

Figure [T]: illustrates the 'Post- X-line Phase' during 
which the particle moves toward the downstream and is 
pitch angle scattered. The rate of energy gain decrease, 
but the particle is energized up to ~0.7mcC^. The ener- 
gization takes place near the edge of the outer diffusion 
region or the outflow exhaust. This is where magnetic 
field magnitude starts to increase. Also, because elec- 
trons are already pre-energized at the X-line, they have 
gyro-radii comparable to the curvature of the magnetic 
field lines of this region. Concurrently, the gradient B 
and curvature B drift motion plays an important role 
for the additional energization during this phase. In the 
later half of the Post-X-line Phase, Particle #1 is re- 
flected by a mirror force at (x, y)~(3.5di, 74d,;). The 
particle again passes through the edge of the diffusion 
region at (a;,?/)~(0, 76di) but the energy increase is very 
small. Finally, it escapes away from the diffusion region. 

We now explore the details of the second rapid ener- 
gization of Particle #1. Figure |8] shows the last three 
of the four phases described in Figure [HI The enlarged 
views of the trajectory are shown in the upper three pan- 
els and the time profile of the particle energy as well as 
the electric field felt by the particle are shown in the 
lower two panels. 

FigurelHK illustrates the 'Pre- Anti-X-line Phase' during 
which the particle approaches toward the merging point 



of two large islands. The particle is reflected twice by the 
mirror force but is not energized. The reversed electric 
field is already large at this phase reaching to ~0.3VaBo 
and the non-magnetized region starts to appear at the 
merging point. 

Figure |Ht> illustrates when the anti-reconnection is tak- 
ing place ('Anti-X-line Phase'). The particle is drawn 
into the anti-reconnection region and is energized signif- 
icantly. The duration of energization, however, is some- 
what longer than the time duration of the first rapid en- 
ergization (ficiAt^6). Note that the anti-reconnection 
is embedded in closed magnetic field lines and the out- 
fiow from the anti- X-line collides with these magnetic 
fields. Moreover, because of the first rapid energization, 
the gyro-radii of the electron is large comparable to the 
scale of the anti-reconnection. The electron motion is 
thus decoupled from the magnetic field lines. The above 
features lead to trapping of the electron within the anti- 
reconnection region and hence significant energy gain. 
Yet, each 'kick' occurs in the inner diffusion region as 
was the case during the first rapid energization. This is 
more clearly shown in the time profiles of the particle 
energy and the x positions (Figure Eli). The oscillating 
feature with the time scale of f2ciAt~l corresponds to 
the gyro-motion of the particle. During the Anti- X-line 
Phase, this oscillatory feature is modified and a contin- 
uous energy increase occurs whenever the cc-position is 
within ±ldi from the X-line. 

Figure (Ht shows the 'Post-Anti- X-line Phase'. Note 
the change of the horizontal axis range. The particle is 
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Fig. 7. — The enlarged view of the trajectory of particle #1 (blue curves) during (a) 62. 0< Qcit <65.7 (Pre-X-line Phase), (b) 65. 7< 
Qci* <68.7 (X-line Phase) and (c) 68. 7< fidi <75.3 (Post-X-line Phase), and the time histories of (d) the particle energy and (e) the 
electric field felt by the particle. The trajectory is superposed on images of (a) the out-of-plane component of the electric field, (b) the 
out-of-plane component of the non-motional electric field and (c) the magnetic field magnitude. The inset shows a blow-up of the X-line 
region of panel b. In each panel, the contour shows the magnetic field lines. The dashed curve in (d) is the same as the solid curve but 
multiplied by 10 to highlight the Pre-X-line Phase. In order to eliminate particle noise, the electric field profiles have been smoothed by a 
box average with the box size of Qci^t^ClS. 
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ejected out of the anti-reconnection region but continues 
to gain energy by the gradient B and curvature B drift 
acceleration. Because of the contracting motion of the 
merged island, the electric field in this region continues 
to increase and reaches 0.4VaBo. As a result, the particle 
energy reaches near 2mcC^. 

In summary. Particle ^1 is first energized at and 
around the X-line generated during the non-linear evolu- 
tion of the tearing mode instability. It is then energized 
at and around the anti- X-line generated by the merging 
of two islands. In both cases, the particle is accelerated 
directly by the reconnection electric field. As a result, 
the energy reaches more than e >mcC^. 

4.2. 0-type Acceleration 

We now show another example of electrons that be- 
haved quite differently from Particle #1 (hereafter Par- 
ticle #2). Figure HI shows the trajectory with the same 
format as Figure [S] but the trajectory is divided into five 
segments instead of four and the background images are 
replaced by the out-of-plane component of the electric 
field instead of the electron current density. 

The first energization occurs at y ^88.5di where an 
X-line develops as a consequence of the first coalescence 
growth phase (Figure [nt-)- There is a small, secondary 
island at y ~86di but it only modulates the particle orbit 
and plays a minor role for the energization. The particle 
travels until it reaches another X-line which is also pro- 
duced by the first coalescence growth (FigurelHb)- There, 
the particle receives two 'kicks' by the reconnection elec- 
tric field. The second kick occurs because the particle 
is reflected by the magnetic fields being piled-up on the 
pre-existing current sheet and comes back to where it can 
again receives energy from the reconnection electric field. 
The two kicks can be identified in Figure [9f at (e/meC^, 
y/di)-(0.5, 45) and ~(1.8, 50). 

What makes Particle #2 different from ^1 is the be- 
haviour during the rest of the trajectory. After the two 
kicks, the particle is trapped within an island and cir- 
culates 6 times (Figure Et) , as can be counted from the 
zigzag orbit of the Az-profile (Figure |9^). Note the fact 
that this island is moving toward the negative y-direction 
because this island and the other island at y ~20di are 
attracting to each other. Thus, the direction of the gra- 
dient B and curvature B drift is positive z at both end 
of the island while the direction of the motional electric 
field is negative and positive at the upper (y ~70di) and 
lower {y ~50di) ends, respectively. As a result, the par- 
ticle gains energy at the upper end but loses energy at 
the lower end. After all, the net energy gain of the parti- 
cle is kept small. This is represented by the foldings and 
overlappings of the curve in Figure [9f at e ^ImeC^. 

The particle eventually reaches the merging point 
where it receives energy at the anti- X-line (Figure |9ji), 
but the anti- X-line has already been well developed so 
that the particle is not trapped in this region. The par- 
ticle soon exits from the merging region and start to cir- 
culate inside the newly created, large island (Figure [St). 
This large island is contracting in this phase so that the 
particle can gain energy every time it passe s through the 
both end of the island (jPrake et al.l 120061 ). There ex- 
ists reversed electric field in between the island edges 
so that the particle can gain energy through the gradi- 
ent B and curvature B drift. As for the displacement 



along the z-direction, the particle moves toward nega- 
tive z-direction only during the period of energization 
by the anti-reconnection. The resultant displacement is 
Az -lOOd,. 

In summary. Particle #2 gains significant amount of 
energy from the motional electric field due to the con- 
tracting motion of the merged island. The final energy 
reaches as high as that of Particle #1, that is e '^2meC^. 

Another important feature of coa lescence is the b ounc- 
ing motion of merge d islands (jTaiima et al. I 119871 : 
I Wan fc Lapmtel 120081 ). The merged island still pos- 
sesses a certain amount of kinetic energy so that a con- 
tracting motion of an island is followed by an expanding 
motion. 

Figure rTOT a-c') shows the evolution of the merged is- 
land. At ilcit=104, the anti-reconnection has stopped 
but the reversed electric field remains because of the 
expanding motion in the x-direction (Figure llOb .) . The 
shape of the merged island is somewhat elongated in the 
y-direction and the contracting motion in the y-direction 
still continues as evident from the negative region. At 
r2cit=112, the merged island has already expanded in the 
x-direction and the contracting motion in the y-direction 
has become weak (Figure [TOb). Notable here is the rip- 
ple structures that surround the island core region, as 
evidently shown by the enhanced, localized E^ regions. 
These structures are due to a turbulent motion within the 
island. The turbulent motion is generated because the 
converging fiows from the top and the bottom are mixing 
with each other. The ripple structures propagate mainly 
in the y-direction away from the merging center. By the 
time r2cit=135, the ripple structure disappears and the 
island has become round. It is very slowly expanding. 

Also shown in Figure fTOl is a sample trajectory of elec- 
trons that pass through a ripple structure (hereafter Par- 
ticle #3). After being energized up to e ~1.3meC^, Par- 
ticle #3 undergoes pitch angle scattering just outside 
an anti-X-line region, i.e. (a:,y)'^(6,38) (Figure [TOk). 
It then travels toward the negative y direction (Figure 
[TUb). During this travel, the particle encounters the 
ripple structure identified at (2;,y)~(6,31). This electric 
field is generated by a local flow toward the positive x di- 
rection and the local magnetic field directed toward the 
negative y-direction. The electron is energized by the 
enhanced E^ in association with the gradient B / curva- 
ture B drift (Figure [TOJl) • A gradient E drift is unlikely 
to be causing the energization because such drift motion 
is perpendicular to E^, the only dominant component at 
this time. 

4.3. Hybrid Type Acceleration 

In a separate paper, we reported an 'island surflng' 
mechanism of electron acceleration during magnetic re- 
connection (jOka et al.l l2010[ ) . This mechanism utilizes 
the secondary magnetic islands that are produced in the 
diffusion region. Inside each island, electrons are trapped 
for a significant period of time so that they are energized 
continuously by the reconnection electric field prevalent 
in the diffusion region. Although a pre-acceleration is 
required to keep electrons trapped within islands, the 
trapped electrons can receive unlimited amount of en- 
ergy as long as the island stays in the diffusion region. 

Note again that while X-Type acceleration takes place 
at any X-line, O- Type acceleration occurs far away from 
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Fig. 8.— The enlarged view of the trajectory of particle #1 during (a) 75. 3< n^t <87.4 (Pre- Anti- X-line Phase), (b) 87. 4< Oci* <93.8 
(Anti-X-line Phase) and (c) 93. 8< fidt <104.0 (Post- Anti-X-line Phase), and the time histories of (d) the particle energy, x positions 
and (e) the electric field felt by the particle. The trajectories in (a-c) are superposed on images of the out-of-plane component of the 
non-motional electric field with the contours showing the magnetic field lines. In order to eliminate particle noise, the electric field profiles 
have been smoothed by a box average with the box size of f2(,i^t~0-13. 



X-lines by taking advantage of a closed field line geom- 
etry. The 'island surfing' mechanism takes place in a 
diffusion region very close to the X-line, but at the same 
time, makes use of the closed field line geometry of a sec- 
ondary magnetic island. As such, this mechanism falls 
into both categories described above and we rather con- 
sider it as an hybrid type. 

The 'island surfing' is a natural consequence of the sec- 
ondary tearing instability of magnetic reconnection, and 
hence not directly related to multi-island coalescence. 
Therefore, we do not discuss its detail any further in 
this paper. 

5. DISCUSSION 

5.1. Acceleration processes 

We performed 2D PIC simulations of multi-island co- 
alescence with no guide field. By following the trajecto- 
ries of 198 energetic electrons, we identified various ener- 
gization mechanisms as have been summarized in Table 
[TJ Note again that these electrons are extracted at the 
peak time of energetic electron fiux with the criteria of 
£ >1.4meC^. The first column shows the categorization 
of each mechanism. The second column shows mech- 
anism nomenclature used throughout this Paper. The 
third column shows the origin of the electric field from 



which particles gain their energies. The fourth column 
shows what we refer to as the 'contribution probabil- 
ity' Pc defined as Pc = WOuacc/ntot where Uacc is the 
number of electrons accelerated by the mechanism and 
ritot is the total number of electrons we analyzed. In the 
present case, ntot=198. For each trajectory of the 198 
electrons, we checked which mechanism the particle had 
been energized through and summed up the number of 
trajectories each mechanism worked on. Because each 
particle experiences more than one different mechanism, 
the total of Pc does not equal to 100%. The fifth col- 
umn shows the largest energy gain by each mechanism 
Ae; . Finally, the last column of Table [1] indicates which 
literature discussed each mechanism. Below, we summa- 
rize each mechanism based on our simulation results and 
describe how we counted riacc for each mechanism. 

The 'surfing' was difficult to identify in our simula- 
tion because of the negligible amount of energy increase 
(Aej/moC^ <0.1). This may be due to the undriven na- 
ture of our simulation that generates polarization elec- 
tric field Ep ^IYaBq. The original study of the 'surf- 
ing' mechanism used a driven magnetic reconnection that 
generates the polarization electric field of Ep ^6Y aBq so 
that the Lorentz force can be balanced by the force from 
Ep. For Table [U we simply did not count the number of 
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Fig. 9.— The trajectory of particle #2 during (a) 46. 7< Qci* <56.0, (b) 56. 0< Q^it <62.0, (c) 62.0< ^ci* <97.3, (d) 97.3< Q^it <99.3 
and (e) 99. 3< f2cit<114.7. Also displayed are (f) the energy and (g) the displacement AZ during 0< ficit <114.7 as functions of Y. The 
background images are the out-of-plane component of the electric field at (a) f2cit=52.7, (b) Qcit=58.7, (c) f2cit=84.0, (d) f!cit=98.0 and 
(e) ncit=104.0. 




-10 10 



-10 10 



-10 10 



idj ' ' ' ' 


III! 


J - 


- - 



1.1 1,2 1,31,4 1.5 1.B 20 24 23 32 



Fig. 10.— The trajectory of particle #3 during (a) 100. 3< r2cit<106.0 and (b) 106. 0< Ccit<114.7. The energy (d) and the displacement 
(e) during these periods are shown in the same panel. The trajectory after f2(.it=114.7 is not shown. The background images are the 
out-of-plane component of the electric field at when the merged island is (a) contracting, f2(.it=104.0, (b) stopped contracting, nt.it=112.0, 
and (c) expanding, f2(,it=135.6. The contour shows the magnetic field lines. 
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the 'surfing' mechanism. 

The 'X-line' mechanism is the classical way of ener- 
gizing electrons, but based on the recently revealed two- 
scale structure, we verified that a significant energization 
takes place at the very center of the diffusion region or 
the 'inner diffusion region' (A£;/meC^ ^0.8). This may 
be a matter of course given the fact that the out-of-plane 
current is localized within the inner diffusion region. The 
number of electrons energized by this mechanism Uacc 
was 117. When counting this mechanism, we made sure 
that the electrons passed through the inner diffusion re- 
gion. 

The 'Anti- X-line' mechanism is the most common ener- 
gization mechanism among the trajectories we analyzed 
(''T-acc=169). It is often intense and and the largest en- 
ergy increment was Ae^/meC^ ~1.8. The energization 
process itself is basically the same as the 'X-line' mech- 
anism but there are three major differences. One is that 
the anti- X-line is created in a driven manner by two sep- 
arate magnetic islands approaching toward each other. 
Since each island is expelled in association with the out- 
flow from an X-line, the anti-reconnection rate reaches 
as large as 1 VaBq. The second major difference is that 
an anti-X-line is bounded by the closed field lines of the 
merged island. The size of the anti-reconnection region 
is relatively small but the electrons can easily be trapped 
at and around the anti-reconnection region. This effect 
leads to multiple number of interaction with the inner 
diffusion region of the anti-reconnection. As a result, the 
'Anti-X-line' mechanism was found to be a significant 
energization mechanism, although the anti-reconnection 
is a transient process that takes place in a small region. 
The third major difference between the X-line and the. 
anti- X-line is that the electric field is reversed at the anti- 
X-line. Because of this, electrons move toward the oppo- 
site direction while being energized. In fact, this feature 
was what we used to count the number of anti- X-line 
mechanism in Table [1] 

The 'grad-B/curv-B drift' mechanism is also a common 
process for electron energization (uacc ~174) but each 
kick is not significant (Ae/moC^ <0.3). Typically, this 
process occurs after a particle exits from either X-line or 
anti- X-line. We identified this mechanism whenever the 
trajectory showed energy increase during a drift motion 
near either X-line or anti-X-line. A drift motion can 
also take place at each end of magnetic island but, in 
such cases, the trajectory can be regarded as a part of 
circulation inside the island and are not counted as the 
'grad-B/curv-B drift' mechanism. The number of the 
'grad-B/curv-B drift' mechanism nacc was 174. 

The 'contracting island' mechanism was identified in 
79 trajectories. Since this mechanism requires energy 
gain at both end of island, the process appears only in 
the later phase of island coalescence. In the earlier phase 
of coalescence, the merging is not complete and either 
of the two approaching island is not contracting. We 
identified this mechanism by the zigzag orbit of the y-e 
plot. The largest energy gain was Ae;/moC^ ~0.6. 

The 'ripple' mechanism was found to be as important 
as the 'contracting island' mechanism, Uacc—^^- During 
the bouncing motion of a newly merged island, a turbu- 
lent flows appear inside the island so that many electrons 
pass through localized electric fields. We counted the 




AZ/di 

Fig. 11. — The displacement of the most energetic electrons along 
the z-direction at different times. The filled circles and the crosses 
are for the electrons in the right and left current sheet, respectively. 
The sign of Az is converted so that electrons of both right and 
left current sheet move toward the same direction as that of the 
primary reconnection electric field. 

number of this mechanism by checking the association of 
an energy increase with the localized electric field within 
a merged island. 

The 'island surfing' mechanism was very rare. Only 4 
electrons were found to be energized by this mechanism. 
This is probably because of the relatively small reconnec- 
tion electric field in the later phase of the simulation run 
when secondary islands are created. Moreover, the size 
of each X-line was limited due to the presence of mag- 
netic islands. If an X-line were able to reach a steady 
state, it would continue to spawn a number of secondary 
islands, which situation would lead to more importance 
of this mechanism. 

In Table [I] we classified each mechanism as either X- 
Type, 0-Type, or a hybrid. An energization mechanism 
belongs to the X- Type if it occurs in the diffusion region 
of magnetic reconnection and/or in the magnetic field 
pile-up region just downstream of the diffusion region. 
The source of energy is the reconnection electric field. 
On the other hand, an energization mechanism belongs 
to the 0-Type if it takes place within the closed geom- 
etry of magnetic field or magnetic islands. The source 
of energy is the kinetic motion of the island. A hybrid 
mechanism refers to an energization within a secondary 
island located in the diffusion region. 

5.2. Displacement along the out-of-plane direction 

The striking feature of electron acceleration at the 
anti-reconnection site is the direction of motion in the 
2-direction being opposite to that of electrons at the pri- 
mary reconnection site. If electrons were magnetized and 
drifti ng, they lose the ir energy by the reversed electric 
field (|Pritchett Il2008f ). However, at the very center of 
the anti-reconnection, electrons are not magnetized and 
not drifting so that they can move parallel to the reversed 
electric field and receive substantial amount of energy. 

Figure [TT] shows the displacement of the most energetic 
electrons along the out-of-plane direction, Az, obtained 
at flcit — 43, 72, 102, and 115. The solid line indicates a 
hypothetical energization by the electric field of O.IVaBq 
which roughly corresponds to the primary reconnection 
electric field measured in the earlier phase of the simula- 
tion. It is also the typical value of steady state magnetic 
reconnection (|Shav et al.,,2007.) . In the earlier phase of 
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TABLE 1 

Electron acceleration mechanisms during multi-island coalescence with no guide field. 



Mechanism Electric field origin Pc,% Aej/moC"^ Reference 



X-Type 


surfing 


reconnection 


N/A 




X-line 


reconnection 


59 




anti-X-line 


anti-reconnection 


85 




grad-B/curv-B drift 


(anti-)reconnection 


88 


O-Type 


contracting island 


island contraction 


40 


ripple 


island ripples 


31 


hybrid 


island surfing 


reconnection 


2 



the simulation run (i7cit=73, green marks), many elec- 
trons follow this solid line because they are energized 
at the primary reconnection regions. They can reach as 
far as 100 di. Later on, electrons are accelerated by the 
anti-reconnection so that they move toward the opposite 
direction (ncit=102, blue marks). By the time of peak 
electron flux (r2cit=115, black marks), the asymmetry of 
the particle distribution remains, but most particles are 
confined within 140 di from the initial positions. 

The fact that many electron positions are largely de- 
viated from the line of energization by the primary re- 
connection, particularly in the later phase of the simula- 
tion run, indicates that these electrons experienced anti- 
reconnection. Some electrons showed unusually large dis- 
placement of Az >140di. All of these particles are ener- 
gized by the 'island surfing' mechanism and so did not 
experience anti-reconnection and direction turnings. 

5.3. Application to the Solar Flares 

Let us now discuss applications of our simulations. We 
used periodic boundary conditions under the assump- 
tion of infinitely long current sheet in the y-direction. 
Therefore, our simulation can probably be best applied 
to the solar flare s which show evidence of very long cur- 
rent sheets (e.g. iSui fc: Holman 1 120031 : iBemporad et al.l 

For the discussion, the energy spectrum in Figure [5] is 
revisited as shown in Figure [Tj] (left). A caveat here is 
that our simulation uses unrealistically high, initial tem- 
perature (eth) due to limited computational resources. 
The temperature of the initial background electrons eth 
= 9.216xl0~3jn^c2 ~ 55MK ~5 keV which is already 
hot in the solar corona. In Figure [12] (left), we re- 
normalized the energy by eth (top axis) and then con- 
verted to the unit of keV under the assumption of a typ- 
ical electron temperature in the solar corona, eth=2MK 
(bottom axis). Moreover, because the typical energy 
range of a space-borne, hard X-ray detector is from a few 
keV to a few hundreds of keV, we only showed the higher 
energy end of the simulated spectrum. For the vertical 
axis, we used the number of electrons F(£) counted in the 
simulation box. F(e) can be representing a population 
n(e) keV~^ cm~^ of electrons in a fully ionized plasma 
of the energy release region. Since this is a simulation, 
the numbers can be scaled arbitrarily. 

As a reference, we reproduced in Figure [12] (right) 
photon spectra from the well-known, solar flare event 
of 23 July 2002 obtained by the RHESSI satellite (e.g. 
iLin et all [2o5l lAsai et al.lf2009D . In the pre-impulsive 
phase, the spectrum showed a power-law-like form that 
extends up to >40 keV from a coronal source. ILin et ahl 
(|2003.) interpreted this spectrum by the combination of 
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an exponential (thermal) and double power-law (non- 
thermal) spectra. In the impulsive phase, the energy 
range below 40 keV is dominated by a single exponential 
form, indicating clearly that electrons are signiflcantly 
heated during this phase. Beyond this 'super-hot' com- 
ponent is the power-law spectrum that extends up to 
>100 keV that originates from the chromospheric foot- 
points of the flare loop. 

It is worth emphasizing here that our simulation is not 
intended to reproduce any particular solar flare event. 
We must keep in mind that our simulation has many 
assumptions and limitations as described below. First, 
the spatial size of the simulation domain is very small 
'-^lOO times the ion inertia length di, that is '-^lO^cm or 
0.14 milli-arcsec. This is much smaller than the typical 
size of the energy release region in solar flares, that is 
^lO^cm or 20 arcsec. Therefore, it is impossible to take 
into account the dynamical evolution of the solar flares. 
Second, the ion to electron mass ratio mi/mo = 25 is 
very small compared to the actual value of 1836. Third, 
our simulation is a two-dimensional model in x and y so 
that spatial variations along the out-of-plane direction z 
are not taken into account. Both reconnection and anti- 
reconnection creates strong current at and around the X- 
lines, and our simulation lacks possible consequences of 
current driven instabilities that would occur in the third 
dimension at the X-lines. Finally, we assumed zero mag- 
nitude of the magnetic field in the out-of-plane direction 
Bz, or the 'guide field'. This assumption is merely for a 
simplicity. In order to better understand the physics of 
the solar flares, we must perform additional simulation 
runs with different guide field magnitude. 

Nevertheless, the fact that our PIC simulation was able 
to produce electrons of up to >30 keV suggests that 
multi-island coalescence is potentially important for the 
understandings of electron acceleration/heating during 
the solar flares. Note that a multi-island coalescence may 
appear in association of developed turbulence that cre- 
ates volume- fllling magnetic island (iTaiima &: Shibata] 
'199/1; iDrake et al.l [20061 : iRetino et al'l 120071 : IBemporad I 
12008( 1. Since turbulence may be generated in any pos- 
sible electron energization site, our simulation may also 
be applied to these sites (Figure [121 bottom). The de- 
tails of the possible electron energiz ation sites as well as 
theories can be found elsewhere (e.g. Asc hwanden |[2002l 
and references therein). It is to be emphasized again that 
we used periodic boundaries but energetic particles that 
reached e >1.4moC^ (or 26 keV when re- normalized) did 
not cross the boundary more than once (e.g. Figure[6]and 
Figure [9|), indicating that electrons do not need to move 
more than lOOdi (or ~10^ cm when re-normalized) in the 
j/-direction to reach the energy. They also do not need 



Electron Acceleration by Multi-Island Coalescence 



13 



to move more than 140di in the z-direction as we showed 
in Figure [11] Note also the total simulation time was 
^130ri~^. If we assume the coronal magnetic field mag- 
nitude to be lOOG, our simulation suggests that electrons 
of up to 30 keV can be created within a spatial extent 
of lO'' cm within a time scale ~ 10/isec. Such a local- 
ized and quick acceleration of electrons is an important 
feature of the multi-island coalescence. 

Our PIC simulation also showed that the total en- 
ergy gained by both thermal and non-thermal electrons 
reaches more than 30% of the magnetic field energy re- 
leased (Figure [3^) . This relatively large number is the 
advantage of having multiple number of X-lines during 
the multi-island coalescence. Our simulation of single 
X-line reconnection with exactly the same parameters 
employed in the presented simulation showed the frac- 
tion of 20% (not shown). Note that X-ray observations 
show that non-thermal electrons alone carry u pwards of 
10-50 % of the released magnetic field energy ([Lin et afl 
120031 and references therein). Therefore, our simulation 
may not fully explain the energy budget during the solar 
flares, but the multi-island coalescence is an attractive 
way to develop future models of the solar flare energy 
budget. 

Figure [T^] suggests that the current PIC simulations 
produce similar spectra as observed in coronal hard X-ray 
sources, but cannot account for the flat electron spectrum 
needed to produce the hard X-ray spectrum observed in 
footpoints. In fact, if we were able to perform PIC sim- 
ulation with a larger simulation domain, say lOOOdiX 
lOOOdi, the third spectral component described in Fig- 
ure [5] may extend to higher energy range to form a clear 
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power-law. Moreover, there would be a variety of dif- 
ferent size of magnetic islands so that the resultant tur- 
bulence leads to a Fermi-like, stochastic acceleration of 
electrons. Therefore, we anticipate our results to be a 
starting point for more sophisticated models of particle 
acceleration during the explosive energy release phenom- 
ena. 

6. CONCLUSION 

We performed 2D PIC simulation to study electron 
acceleration during multi-island coalescence. By ana- 
lyzing the trajectories of the most energetic electrons, 
we found a variety of different acceleration mechanisms 
such as the contracting island mechanism, ripple mech- 
anism, island surfing mechanism, etc. However, a sta- 
tistical study showed that the most important process 
is the energization process that takes place at the anti- 
reconnection region. Based on the maximum energy of 
electrons attained in the simulation, we pointed out that 
the multi-island coalescence may play an important role 
in producing energetic electrons during the solar flares. 
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Fig. 12. — (top left) Energy spectrum of electrons obtained by the PIC simulation. It is exactly the same as the one in Figure [5] but 
only shows the higher energy end. The background temperature of 2MK is assumed to derive the electron energy. See texts for other 
assumptions and limitations of the simulation model, (top right) Energy spectra of X-ray photons obtained by the RHESSI satellite on 
2002 July 23. The upper and the lower sp ectra are obtained at 0021:42 UT and 0028:10 UT, respectively. These are the reproduction of 
the energy spectra reported in Figure 2 of I Lin et al.l pOOSl which paper contains the details of the event, (bottom) Schematic illustration 
of possible electron energization sites as well as typical coronal and footpoint sources of hard X-rays. 



